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Abstract 



New models have recently been proposed in which a second Higgs doublet couples only to the lepton doublets 
and right-handed neutrinos, yielding Dirac neutrino masses. The vacuum value of this second "nu-Higgs" 
doublet is made very small by means of a very softly-broken Z2 or U(l) symmetry. The latter is technically 
natural and avoids fine-tuning and very light scalars. We consider a supersymmetric version of this model, 
in which two additional doublets are added to the MSSM. If kinematically allowed, the decay of the heavy 



MSSM scalar into charged nu-Higgs scalars will yield dilepton events which can be separated from the 



W-pair background. In addition, the nu-Higgsinos can lead to very dramatic tetralepton, pentalepton and 



\ hcxalepton events which have negligible background and can be detected at the LHC and the Tevatron. 
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I. INTRODUCTION 



One of the simplest extensions of the Standard Model, and one of the most studied, is the 
addition of an additional Higgs doublet There are many versions of the Two-Higgs Doublet 
Model (2HDM) in the literature, which differ in the couplings of the Higgs doublets to fermions. 
The most familiar are Model I, in which one doublet couples to all fermions and the other couples 
to no fermions, and Model II, in which one doublet couples to the Q = 2/3 quarks and the other 
couples to the Q = —1/3 quarks and leptons. The latter is a feature of supersymmetric models 
Q]. These couplings are generally restricted by a Z2 symmetry. Another model, Model III, has no 
such discrete symmetry but then also has tree- level flavor-changing neutral currents Q]. 

Another version of the 2HDM has one doublet coupling to all of the quarks and the second 
doublet coupling to the leptons. Although the basic structure of this model was proposed long ago 
[4], it has received a resurgence of interest [jj due to the existence of non-zero neutrino masses. 
A modified version, in which one doublet couples to all of the quarks and charged leptons, and 
the second doublet couples only to the left-handed lepton doublet and the right-handed neutrino 
was proposed by Gabriel and Nandi [g]. The model will allow for Dirac neutrino masses, which 
are small due to a very small (less than an eV) vacuum expectation value for the second doublet 
(as opposed to the one-doublet case in which they are small due to very small Yukawa couplings). 
The model has a Z2 symmetry in which the second doublet and the right-handed neutrinos are 
odd and all other fields are even. 

The model of Gabriel and Nandi [6] has the feature that it contains a very light scalar which 
causes problems with standard cosmology. This feature remains even if the Z2 symmetry is pro- 
moted to a U(l) (which can eliminate Majorana mass terms for the right-handed neutrinos). Very 
recently, Davidson and Logan (DL) proposed [7] enforcing the coupling structure with a global 
U(l), but breaking the U(l) explicitly, through a dimension-2 soft term in the Higgs potential. 
This avoids any Goldstone bosons and other light scalars, and only requires that the soft term 
have a magnitude of approximately an MeV 2 . Since this term is the only U(l) breaking term, 
the smallness of its size is technically natural. The charged Higgs boson in the model has very 
distinctive signatures, decaying into a left-handed charged lepton and a neutrino, with branching 
ratios determined by the Pontecorvo, Maki, Nakagawa, Sakata (PMNS) neutrino mixing matrix. 

In this paper, we consider the supersymmetric version of the Davidson-Logan model. We extend 
the MSSM by adding two additional Higgs doublets of opposite hypercharge, which have opposite 
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quantum numbers under the global U(l). The right-handed neutrinos are also charged under the 
U(l) and all other fields are neutral. We first show that, through D-terms, the heavy neutral 
scalar of the MSSM will decay (if kinematically allowed) into the lightest of the [/(l)-charged 
Higgs, leading to distinctive signatures. We then look at the supersymmetric partners of the 
J7(l)-charged Higgs, and find some remarkable signatures at the LHC and the Tevatron, including 
dramatic multilepton events which have little to no background and yet accessible cross sections. 

II. THE MODEL 

To ensure anomaly cancellation, we must add an even number of doublets to the MSSM. Re- 
ferring to the MSSM doublets as Hi and H2, we add two new doublets, H3 and H4 and three 
right-handed neutrino superfields, JVj. A global U(l) symmetry is imposed in which H3, H4 and 
Ni have charges —1, +1 and —1 respectively, while all other fields are uncharged. As usual, R-parity 
will be conserved in the model. 

The superpotential is 

W = W + y'NLH 4 + fi 34 H 3 H 4 (III) 

where Wo is the MSSM superpotential and L is the left-handed lepton doublet (generation indices 
are understood). Since the new Higgs fields only couple to right-handed neutrinos and leptons, we 
refer to them as "nu-Higgs." Due to the global U(l) symmetry, there is no mixing with the MSSM 
Higgs bosons at this stage. We also assume that the mass-squared parameters of the nu-Higgs 
fields are sufficiently large and positive that the fields do not acquire vacuum expectation values 
(vevs). DL show that there is a cosmological lower bound of 1/30 on y'; our results will not be 
sensitive to y' as long as it is not extremely small (so that decays occur at the vertex). 

Following Davidson and Logan, the global U(l) is broken explicitly by soft dimension-2 terms 
in the potential Hi^HiH^ + ^33-^2^3- These terms will result in vevs for the nu-Higgs fields which 
have the seesaw form V3 = ^2^2 /Ma and V4 = (J-i^vi/Ma, where Ma is generally the weak scale. 
Since the vevs of H3 and H4 cannot be too dissimilar (due to the [/(l)-conserving /i 34 coupling), 
we must have V3 ~ V4 ~ eV in order to give the correct neutrino masses without very small Yukawa 
couplings. This implies that and /X23 are both of the order of an MeV. 

Since these soft terms are the only source of U(l) breaking, their small size is technically 
natural. [8] One can easily imagine getting such small terms through a Froggatt-Nielsen mechanism. 
For example, introducing a flavon field F with a U(l) charge of 1/3 would allow one to include 
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an effective U(l) invariant operator j^F^/j^^Hs, where M is the flavon scale and \x is the weak 
scale. Then for (F)/M ~ 0.05, one obtains an effective term of the correct size. We will not pursue 
this issue further, but just note that it may not be too difficult to explain the size of these small 
terms. 

This model does seem rather ad hoc, and one can ask about a possible ultraviolet completion, 
possibly in the context of a grand unified theory. Of course, this question can be asked about any 
leptophilic model, whether super symmetric or not. As noted earlier, there are many discussions of 
leptophilic models HQ 

in which the charged leptons (and the left-handed neutrino) have a different 
sign under the Z2 symmetry than the quarks. It is difficult to see how to reconcile these models with 
grand unification, since the charged leptons and quarks are generally in the same representation. 
In this model, however, the charged leptons and quarks have the same sign under the Z2, and thus 
it is easy to incorporate into a GUT. For example, in SU(5), one could simply promote the H fields 
to 5-plets, and the N fields to singlets, and impose the same discrete symmetry. Whether one can 
find a model which can also give Froggatt-Neilsen type terms discussed in the previous paragraph 
is unclear — it would possibly be more contrived than the model we have presented. Without such 
a detailed model, discussion of sparticle mass spectra would be premature. Nonetheless, this model 
would appear to be easier to embed in a GUT than other leptophilic models. 

The nu-Higgs spectrum consists of two scalars, two pseudoscalars, two pairs of charged scalars, 
two neutral nu-Higgsinos and a pair of charged nu-Higgsinos. What are their masses? In the scalar 
sector, DL had several unknown quartic couplings, and thus the relative masses of the scalars was 
arbitrary; they had to consider cases in which the charged scalar was either heavier than or lighter 
than the neutral scalars, which affected the phenomenology. In this case, however, the quartic 
terms are completely determined by gauge couplings. Ignoring the very small corrections due to 
U(l) breaking, the masses were calculated in Ref. [{|. They found that the neutral scalar mass 
matrix depended on unknown parameters, but the pseudoscalar mass matrix was identical to the 
scalar, and the charged scalar mass matrix only differed slightly. It was shown that the lightest 
charged scalar was a few GeV heavier than the degenerate scalar and pseudoscalar. Since the 
masses are similar, decays of the nu-Higgs particles into each other will be phase-space forbidden 
or suppressed, and the decays into leptons will predominate. As a result, we will focus on the 
lighter of the states, and take the scalar, pseudoscalar and charged scalar to be degenerate in mass 
(for phenomenological purposes, keeping in mind that the charged scalar is slightly heavier). They 
will be referred to as x°>XA and respectively. The nu-Higgsino masses depend only on ^34 
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and are completely degenerate at tree level. Note that the neutral nu-Higgsino mass matrix is 
completely off-diagonal, giving a mixing angle of 45 degrees when rotating to mass eigenstates. 
Thus, both of the neutral nu-Higgsinos as well as the charged nu-Higginos are degenerate in mass. 
For notational simplicity, we refer to them as x° and X s1 . Although there are two neutral states, 
the branching ratios will be the same for each, and we will account for the factors of two in the 
production cross section. 

Focusing on the lighter of the nu-Higgs and on the nu-Higgsinos, we can look at their decays. 
The decays of the neutral nu-Higgs, X°>XA are into neutrinos, and are thus unobservable. The 
charged nu-Higgs, x , wm decay into all nine combinations liUj. Davidson and Logan show that 
the decay rate into l^u is proportional to £N m1.\Uu\ 2 , where U is the PMNS mixing matrix. For 
a normal neutrino mass hierarchy, with the lightest neutrino having a mass below 10 -3 eV, the 
decays will be into \xv and tv with branching ratios between 40 and 60 percent. For an inverted 
hierarchy, the rate into ev is about 50 percent, with the balance shared equally between [iv and 
tv. A discussion of the PMNS mixing matrix parameters used is in Appendix [A] The actual width 
is narrow, but the decay still occurs at the vertex. 

The nu-Higgsinos decay as 



with equal branching ratios in the limit that the slepton and sneutrino have equal mass. Note that 
because the nu-Higgsino mass eigenstates are made up of equal parts of H% and H4, and the former 
does not have any two body decays, both mass eigenstates will decay similarly. 

Finally, the right-handed sneutrino, vr can only decay through a virtual nu-Higgs (plus a left- 
handed sneutrino or slepton) or a virtual nu-Higgsino (plus a left-handed neutrino or lepton), and 
has no two body decays. As a result, it can have visible decays, and thus the x° decays above 
could be observable through both decay chains. 

We now turn to the detailed phenomenology of these nu-Higgs particles at the LHC and the 
Tevatron. 

III. PHENOMENOLOGY OF THE NU-HIGGS SCALARS 

In the non-super symmetric version of the model, the phenomenology of the nu-Higgs scalars 
was discussed by Davidson and Logan [7|. As noted above, since supersymmetry did not restrict 




(II.2) 
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the masses in their model, they needed to consider a range of possibilities. In the super symmetric 
version, the lightest neutral nu-Higgs is only slightly lighter than the lightest charged nu-Higgs. 
As a result, it will only decay into neutrinos and would thus be unobservablefisl] . The charged 
nu-Higgs will decay into a charged lepton and a neutrino; summing over the unobserved neutrino 
leads to the charged lepton being approximately 50% fx, 50% r (50% e, 25% fj,, 25% r) for the 
normal (inverted) hierarchy. 

This was discussed by Davidson and Logan, and the primary production mechanism was through 
a Drell-Yan photon or Z. In the supersymmetric version, a new production mechanism opens up, 
and can be substantially larger. Because of the 577(2) and U(l) D-terms, there is an interaction 
(h\Hi — ff| -^2) (-£^3 -£^3 — H\H/l), and when H\ and H2 acquire vev's, this will lead to a three point 
interaction between an MSSM Higgs and two nu-Higgs scalars. The lightest MSSM Higgs is too 
light to decay into two charged nu-Higgs, but the heavier one, H, is not. One can thus produce 
the H of the MSSM via gluon fusion [101 ]. and it will then decay into a pair of charged nu-Higgs 
scalars. Each of those will decay into a charged lepton and a neutrino. 

For the production of the heavy MSSM Higgs through gluon fusion we use tan f3 = 3. In order 
to calculate the pair production of the charged nu-Higgs bosons we first calculate the branching 
fraction of the heavy MSSM Higgs to a pair of charged nu-Higgs bosons. From the D-terms in the 
potential we can derive the coupling of the heavy MSSM Higgs to the charged nu-Higgs. 



where a and j3 are the standard MSSM parameters and 7 is the mixing angle that diagonalizes the 
neutral nu-Higgs mass-squared matrix. The decay width of the MSSM Higgs to charged nu-Higgs 
bosons then can be written as 



For example using an MSSM Higgs of 400 GeV and an charged nu-Higgs mass of 100 GeV we 
calculate the production cross section of a pair of charged nu-higgs bosons to be 14 fb assuming 
cos(27) = 1. Since both of the charged nu-Higgs bosons will each decay into a charged lepton and 
neutrino, the collider signature will be two leptons (which generally can be different) and missing 
energy. 

There are two backgrounds that we consider. The primary background is VF-pair production 
with both W bosons decaying leptonically. We also consider tt pair production as another possible 
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background. Due to its large production cross section (~ 0(100 pb)) it is possible for it events to 
have both tops decay semi-leptonically and both 6-quarks be missed by the detector because of too 
large rapidity or too small transverse momentum (i.e. rjb > 2.0 and pT,b < 20 GeV, respectively). 
We generated events for the signal and both of the backgrounds using the MadGraph/MadEvent 



program package ll| using the CTEQ6L parton distribution functions [12]. To simulate these 
events in a collider we apply the following acceptance and isolation cuts on the two final state 
leptons: 

PTMn = 20 GeV, (III.5) 
\ri e \ < 2.4, < 2.1, (III.6) 
AR = VAr/ 2 + A<p 2 > 0.4, (III.7) 

where rj and (j) are the pseudorapidity and azimuthal angle of a final state lepton respectively. 
We also apply a 90% tagging efficiency for the leptons (electrons and muons). For the final state 
leptons we do not include taus because of the small 40% tau tagging efficiency and a relatively 
large light jet mis-tagging rate ~ 0(1%) \}\ - 

To extract the signal from our background we first require that there are two final state leptons 
that pass the acceptance and isolation cuts. This naturally reduces the number of ti events. 
Because our final state particles decay from a heavy Higgs boson, we expect that there should be 
a large amount of missing transverse energy from the two energetic neutrinos as well as a large 
total transverse mass from the high px leptons. We exploit this by requiring the missing transverse 
energy be greater than 90 GeV and the total transverse mass be greater than 250 GeV. Applying 
these collider cuts using the MadAnalysis program package ll| we obtain the results given in 
Table [H With 100 fb _1 of integrated luminosity we obtain a statistical significance near 3a for 
both values of the charged nu-Higgs mass and a heavy MSSM Higgs mass of 400 GeV. For an MSSM 
Higgs of 300 GeV and charged nu-Higgs mass of 100 GeV we obtained a statistical significance over 
5<7 due to the larger production cross section of the MSSM Higgs and larger branching fraction 
to a charged nu-Higgs pair. Note we did not assume any particular hierarchical structure of the 
neutrino masses. If we knew more specifically what the flavor distribution of the leptons in the 
charged nu-Higgs decay, it may be easier to extract the signal from backgrounds. 
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Cut 


w+w~ 


tt 


M H 

M x 


= 400 GeV 
= 100 GeV 




M H 

M x 


= 400 GeV 
= 150 GeV 




M H = 300 GeV 
M x = 100 GeV 


6' 


VS+B 


VS+B 


VS+B 


tagging: 21, jets 


113000 


610 




420 


1.2 ±0.06 




300 


0.9 ±0.05 


2200 


6.5 ±0.1 


E T ,miss > 90 GeV 


1700 


160 




120 


2.7 ±0.2 




150 


3.3 ±0.3 


400 


8.4 ±0.4 


M T r > 250 GeV 


1300 


160 




110 


2.8 ±0.3 




130 


3.3 ±0.3 


210 


5.1 ±0.3 



TABLE I: Collider cuts used to extract the two lepton, two neutrino signal from x + X~ production from the 
backgrounds considered. With 100 fb _1 of integrated luminosity we obtain a statistical significance over 5<r 
for the Mh = 300 GeV case and near 3cr for Mh — 400 GeV and both values of M x . It is interesting to 
note that for the charged nu-Higgs of 150 GeV that even though the production cross section is smaller, 
the larger nu-Higgs yields events with larger missing energies and total transverse mass which allows more 
events to pass our collider cuts. 



IV. PHENOMENOLOGY OF THE NU-HIGGSINOS 
A. Production 

The nu-Higgsinos cannot be produced from decays of an MSSM Higgs, thus the primary pro- 
duction mechanism is through a Drell-Yan process. The charged nu-Higgsino pair is produced 
through a Drell-Yan process involving a 7 or Z boson while the neutral nu-Higgsino pair can only 
be produced through a Z. Note that the Z will only couple to two different nu-Higgsinos, but since 
the two neutral nu-Higgsinos are degenerate this will effectively be a pair production mechanism. 
Associated production of a charged and neutral nu-Higgsino can occur through a Drell-Yan process 
involving the W bosons. 

The LHC and Tevatron production cross sections of the above channels are calculated using the 



MadGraph/MadEvent software package [ll| with the CTEQ6L parton distribution functions [12]. 
Using center of mass energies of 14 and 7 TeV for the LHC and 2 TeV for the Tevatron we calculated 
the Higgsino pair production cross sections for Higgsino masses within the range of 100 — 500 GeV. 
The results are given in Fig [U For the LHC, the x + ~x° P air nas the largest production cross 
section. Note that the x + 'X° cross section is larger than the x~ /x° production due to the valence 
up-quarks used in the former process as opposed to the valence down-quark in the latter. We 
also note the larger charged nu-Higgsino pair production cross section compared to the neutral 
nu-Higgsino production from the Drell-Yan process involving a photon. 

B. Decays 

We first discuss the detailed decays of the neutral and then the charged nu-Higgsinos. For 
simplicity, we will assume that the left and right handed sneutrinos have similar mass. The x° 
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LHC - 14 TeV 
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i.i ' 1 ' 1 ' 1 ' 1 
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nu-Higgsino Mass (GeV) 
(a) 




200 300 400 

nu-Higgsino Mass (GeV) 



(b) 



Tevatron - 2 TeV 






+/- 


ppbar 


->x X 


ppbar 


-> X + X 


■ ■ ■ ppbar 


u 

->xx 



300 

nu-Higgsino Mass (GeV) 



(c) 



FIG. 1: Production cross sections of various pairs of charged and neutral nu-Higgsinos (in fb) at the LHC 
with center of mass energies of (a) 14 TeV and (b) 7 TeV as well as the Tevatron with a center of mass 
energy of (c) 2 TeV. 
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decays into Branch A: vr + T>l or Branch B: i>l + vr . Since the mass of the neutrinos arises from 
the same Yukawa coupling term as this decay, the coupling is flavor diagonal and proportional to 
the mass. 

Let us first consider Branch A. The left-handed sneutrino decays are well-studied in the MSSM. 
Here, we must choose a point in the mSUGRA parameter space. We will choose Snowmass point 
SPS-la, which is one of the most studied to date. Using the results of Freitas et al we find that, 
for the normal hierarchy, the left-handed sneutrino decays into Bvl 88% of the time, either vlt + t~ 
or Bult + t~ 7% of the time, vu^ 1 ^ 5% of the time, and has a negligible (0(0.2%)) decay into 
muon or electron pairs. For the inverted hierarchy, the left-handed sneutrino decays into Bvl 88% 
of the time, either vlt + t~ or Bult + t~ 5% of the time, ^r^e^ 5% of the time, i/lt ± /x =f 2% of 
the time, and also has a negligible (0(0.2%)) decay into muon or electron pairs. Since only 12% 
of the decays are visible in either hierarchy, the branch into the left-handed sneutrino is not as 
promising as for the other branch. 

Let us now consider Branch B: The right-handed sneutrino decays have not been previously 
studied. Since there are no gauge interactions, and negligible mixing with MSSM fields, the only 
decays of the right-handed sneutrino are three-body decays via a virtual nu-Higgsino. One can 
have Branch Bl: vlX°* or Branch B2: lx + *, where the * superscript indicates a virtual state. 
Branch Bl will give exactly the same flavor structure as Branch A, and thus we use those results. 
For Branch B2, however, the final state will always be visible, leading to the most interesting 
signatures. 

Consider the normal hierarchy. The right handed sneutrino has a specific flavor, and thus the 
flavor of the lepton in the decay leads to a muon 50% of the time and a tau 50% of the time (see 
Appendix A for a discussion). As before, the x + * wm decay with the same 50-50 split. Thus, one 
expects Branch B2 to give fxjl, ht,tJ1 and rf each 25% of the time. For the inverted hierarchy, 
the flavor of the lepton in the decay of the right handed electron sneutrino is always that of an 
electron, and the decay of the right handed muon or tau sneutrino is into a muon 50% of the time, 
and tau 50% of the time. The x + * a ^ so decays with the 50-25-25 split in flavor so that branch B2 
in the inverse hierarchy yields ee, ejl, er, fie, fjjl, frr, re, tJI, and rf with respective probabilities 
of 25%, 12.5%, 12.5%, 12.5%, 6.25% J5.25%, 12.5%, 6.25%, and 6.25%. Finally, the slepton decays 
are given in Table 1 of Freitas et al. [LJ] . 

Putting all of these together, we find the decays of the x° given in Table [TT] and Table IIIH We 
have only included decays whose branching ratios exceed 0.5%. Of course, each of the r leptons 
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will decay 17% of the time into electrons and 17% of the time into muons. 
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TABLE II: Decay branching fractions for the nu-Higgsinos assuming a normal hierarchy. For the charged 
mi-Higgsino, the lepton in parenthesis has the upper sign for y + and the lower sign for 



We now turn to the decays of the y + - Here our work is already done from the discussion of 
Branch B2 above. The resulting branching ratios are given in Tables llTl and ITTT1 Note that one can 
have single lepton, tri-lepton or penta-lepton decays. 



















Signature 


Branching Fraction 


Signature 


Branching Fraction 


Invisible 


66% 




28% 




11% 


(e±) 


23% 


T + r - 


6.7% 




11% 




5.1% 


(e ± )r+r- 


5.9% 


e + e~ 


3.1% 


(r ± )e+e- 


4.1% 




3.0% 


( M ± )r+r- 


3.6% 




1.1% 


(e ± )e+e- 


3.3% 




1.1% 


(r^)(e ± e ± ) 


3.2% 


T + T - e ± T T 


1.0% 


(r ± )r+r- 


2.9% 




.91% 




2.1% 






(M ± )e+e- 


1.7% 






( e T)( T ± T ±) 


1.1% 






(e ± )e + e t+t 


1.1% 






(r±)/i+/i- 


0.84% 






(^ )(e ± e ±) 


0.75% 






( e ±) M + At - 


0.69% 






(// ± )e + e t + t 
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TABLE III: Decay branching fractions of nu-Higgsinos assuming an inverted hierarchy. For the charged 
nu-Higgsino, the lepton in parenthesis has the upper sign for y + and the lower sign for y~. 
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C. Signatures 



The high multiplicity of charged leptons in the nu-Higgsino decays gives some remarkable signa- 
tures. Prom pair production of the x + X~ i one can have very dramatic hexalepton and tetralepton 
events. Associated production of X s1 with x° leads to pentalepton events, and pair production of 
~o~o i eac i s ^ tetralepton events. The production cross sections were given earlier. 

1. Hexalepton events 

Perhaps the most dramatic events are hexalepton events which arise from the decay of a charged 
nu-Higgsino pair. From Tables HT1 and iHlj one can see that the x has three charged leptons in the 
decay 31 — 32% of the time, for either hierarchy, and thus the pair will yield six leptons roughly 
10% of the time. Note that the production rate at the 14 TeV LHC for a light X s1 P arr can De as 
large as 2600 fb, leading to an enormous hexalepton rate of over 200 fb. 

One will get a distribution of lepton flavors depending on the hierarchy. How robust are these 
results? If the x* 1 is lighter than the slepton or sneutrino, the flavor structure will not change- 
it simply means that the slepton or sneutrino is virtual, thus these results won't change. We 
did choose a specific point in the mSUGRA parameter space, and that will change the flavor 
distribution, and thus these precise percentages should be taken cum grano salis. Nonetheless, one 
does expect a huge production rate for hexalepton events if the x^ 1S n °t too heavy. 

In Table [IV] we have summarized the signatures as follows. In the first row, for both hierarchies 
and X s1 masses of 100, 200, 300 GeV, we give the total production rate for hexalepton events. The 
tau's in the decays will either decay leptonically or hadronically. If they decay hadronically, roughly 
40% will be identified. If they decay leptonically, they will be indistinguishable from a muon or an 
electron. In the remaining rows, we have listed the possible signatures, depending on the number 
of hadronic tau-identifications, weighted by the 40% factor. Since the electrons and muons are 
detected with virtually 100% efficiency, all of the decays listed in Table UVl are detectable. 

One sees that each of the various signatures can be well within reach of the LHC, and some can 
have production cross sections of many tens of femtobarns. We assumed here that \fs = 14 TeV. 
If it is 7 TeV, then the above figures show that the production rates are lower by a factor of 
roughly 3. At the Tevatron, one can scale the results as shown in the production cross section 
figures. One might be able to extend the sensitivity by consider the explicit charges. For example, 
one combination in the normal hierarchy case would be n + fi + /j> + fi~ t~ t~ which might have lower 
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TABLE IV: Cross sections, in femtobarns, for hexalepton signatures for various values of the charged nu- 
Higgsino mass and for the two different neutrino mass hierarchies. The first line gives the total event rate 
into six leptons. For the other lines, we include the fact that leptonic r decays lead to additional /z's and 
e's and that hadronic r decays are detected with a 40% efficiency. In the table t refers to either muons or 
electrons, and t refers to an identified hadronically decaying tau. The quoted cross sections are assuming 
a center of mass energy of 14 TeV at the LHC. A center of mass energy of 7 TeV reduces the cross section 
roughly by a factor of 3. 



backgrounds compared with three charge pairs. Given the number of combinations, a detailed 
analysis of these possibilities would be premature. Note that the total rate of identified hexalepton 
events, for a x ± mass of 100 GeV is approximately 50 femtobarns at the LHC for y/s = 14 TeV and 
approximately 16 femtobarns for = 7 TeV. At the Tevatron, the rate would be roughly 4 fb, 
leading to a couple of dozen events for the current integrated luminosity. Note that the production 
rate at the Tevatron will drop more quickly with the mass than the LHC. 

Suppose hexalepton events are not detected? It could simply be that the charged nu-Higgsinos 
are too heavy, and thus failure to detect their decays would effectively place a lower limit on their 
mass. Could a 100 GeV charged nu-Higgsino evade detection? The choice of hierarchy and the 
choice of the particular mSUGRA point will affect the flavor distribution (which can thus affect 
the tau content of the events and thus the efficiency of detection), but those will not substantially 
affect the size of the signal. If the sneutrinos are heavy, this will also not affect the decays, unless 
they are so heavy that charged nu-Higgsino decays into a nu-Higgs plus the LSP can dominate, 
leading to only two leptons in the decay. Perhaps the simplest way to reduce the signal is to have 
the right-handed sneutrino be substantially heavier than the left-handed sneutrino, since the latter 
decays are much more likely to be invisible. Nonetheless, for most of parameter space, the charged 
nu-Higgsino mass is the primary factor in the signal rate. Note also that, independent of any of 
these factors, failure to detect these events at the Tevatron would imply a lower bound on the 
necessary luminosity to detect them at the LHC, as one can read off from the production cross 
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sections. 

2. Tetralepton and pentalepton events 

X pair production can also lead to tetra-lepton events, as well as x° P a i r production. In 
Table |V] we have listed the total production rate, as well as the various signatures, for these 
events. Finally, associated production from a W will lead to pentalepton events, and these are also 
shown in Table IVII 

The cross sections are larger than for hexalepton events, however backgrounds will also be 
larger. For example, a tetralepton signature of fi + fi~e + e~ would have backgrounds from real or 
virtual 7's and Z's. One would expect the pentalepton signature to have smaller backgrounds, and 
yet the cross sections are still much larger. They are much more sensitive to the flavor distribution, 
as one can see in the table. It may very well be that the pentalepton signatures will be the most 
sensitive. The cross sections are somewhat lower than for tetraleptons, but one would expect the 
backgrounds to be substantially lower. The hexaleptons are more dramatic, with even smaller 
backgrounds, but the cross sections are substantially higher for the pentalepton events. 

How sensitive is the Tevatron? Consider a nu-Higgsino mass of 200 GeV. The cross section for 
pentaleptons at the LHC for ^/s = 14 TeV is approximately 50 fb, as seen in Table IVT1 Comparing 
cross sections between Figure and [lb, one sees that the cross section at the Tevatron is roughly 
a factor of 50 smaller, giving a cross section of 1 fb. Thus, with 8 fb -1 integrated luminosity, one 
would expect 8 events, iff all r's could be identified. Obviously, they can't, and Table IVII shows 
that events with five non-tau's occur at a rate of roughly 20% of the total. Of course, we have not 
done a detailed analysis. Some hadronic r's would be identified, and some of the leptons might not 
be identified (some would not have enough transverse momentum, for example). Although we are 
not aware of any backgrounds that large, a full simulation should be carried out. It does appear 
that a lower bound of approximately 200 GeV on the nu-Higgsino masses could be obtained at the 
Tevatron. Certainly, at a mass of 100 GeV, the event rate is a factor of 10-20 larger and detection 
should be straightforward. 

V. CONCLUSIONS 

The discovery of nonzero neutrino masses has led to a new type of two Higgs doublet model 
in which a second Higgs doublet couples only to the lepton doublets and right-handed neutrinos, 
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TABLE V: Cross sections, in femtobarns, for tetralepton signatures for various values of the charged nu- 
Higgsino mass and for the two different neutrino mass hierarchies. The first line gives the total event rate 
into four leptons. For the other lines, we include the fact that leptonic r decays lead to additional /u's and 
e's and that hadronic r decays are detected with a 40% efficiency. In the table i refers to either muons or 
electrons, and t refers to an identified hadronically decaying tau. The quoted cross sections are assuming 
a center of mass energy of 14 TeV at the LHC. A center of mass energy of 7 TeV reduces the cross section 
roughly by a factor of 3. 





Normal Hierarchy 


Inverted Hierarchy 




100 GcV 


200 GcV 


300 GcV 


100 GcV 


200 GcV 


300 GcV 


5 Lcptons 


770 


64 


13 


600 


50 


10 


Or U 


29 


2.4 


0.50 


110 


9.0 


1.8 


It U 


49 


4.1 


0.84 


100 


8.3 


1.7 


2t U 


36 


3.0 


0.61 


41 


3.4 


0.69 


3t 21 


12 


0.97 


0.20 


8.9 


0.74 


0.15 


At U 


1.8 


0.15 


0.031 


0.76 


0.063 


0.013 


5t m 


0.11 


0.0090 


0.0020 


0.015 


0.0010 


0.000 



TABLE VI: Cross sections, in femtobarns, for pentalepton signatures for various values of the charged nu- 
Higgsino mass and for the two different neutrino mass hierarchies. The first line gives the total event rate 
into five leptons. For the other lines, we include the fact that leptonic r decays lead to additional /x's and 
e's and that hadronic r decays are detected with a 40% efficiency. In the table I refers to either muons or 
electrons, and r refers to an identified hadronically decaying tau. The quoted cross sections are assuming 
a center of mass energy of 14 TeV at the LHC. A center of mass energy of 7 TeV reduces the cross section 
roughly by a factor of 3. 



leading to Dirac neutrino masses. We consider a super symmetrized version of this model, which 
thus contains four Higgs doublets. A global symmetry which is only very weakly broken prevents 
substantial mixing between the additional doublets, called nu-Higgs multiplets, and the MSSM 
Higgs sector. 

We study two aspects of the phenomenology of these new states. For the scalar nu-Higgs 
fields, a new production mechanism leads to the possibility of a 5<r detection with 100 fb _1 of 
integrated luminosity. The second aspect contains extremely exciting possibilities. For the nu- 
Higgsino states, we find remarkable phenomenological signatures. The most dramatic of these 
signatures are hexalepton events, which contain six charged leptons and missing energy and which 
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are produced with a cross section of up to 250 fb for a nu-Higgsino mass of 100 GeV. Many of 
these leptons are tau's, and when we fold in the 40% detection efficiency, we find roughly 50, 5, 1 
fb cross sections for events in which all six leptons can be identified for nu-Higgsino masses of 
100, 200, 300 GeV respectively. For \fs = 7 TeV, these numbers are lower by a factor of three, 
and for the Tevatron are lower by another factor of four. Thus, this model gives the exciting 
possibility of substantial multi-lepton events which can be detected at the LHC and may, for lower 
nu-Higgsino masses, be detectable at the Tevatron. 
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Appendix A: PMNS Mixing Matrix Parameters 

The decay rate of \ + ( as wen as the rate for supersymmetric versions) is proportional to 

the square of the ^-element of the neutrino mass matrix in the flavor eigenstate basis. It, in turn, 
is related to the physical neutrino masses through the PMNS matrix Uij by mf? = Uj k m Uk[ Uij , 
where mf? is the diagonalized neutrino mass matrix. This gives 

nx + -> t+vj) oc m\. = (u ik m° ki Ul) 2 . (A.l) 

The diagonalizing matrix and diagonalized mass matrix are given by (in the limit where #13 = 0) 



U, 



cos(6>i 2 ) sin(0i 2 ) 

- sin(0i 2 ) cos(6> 23 ) cos(6>i 2 ) cos(0 23 ) sin((9 2 3) 

^ sin(6>i 2 ) sin(0 23 ) - cos(6>i 3 ) sin(0 23 ) cos(0 23 ) J 
( 



(A.2) 



and vnP,_ . 



m 



m 



x ' 





V 







m 



(A-3) 



v?, 



I 



Although there is a range of allowed values for #i 2 , 6*13, and 623, we will choose them to be given 
by 34°, 0°, and 45°, respectively. We will also assume that the lightest neutrino has negligible 
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mass (thus not considering the fully degenerate case). The reason for these assumption is that, 
in determining branching ratios, we are already choosing a specific point in mSUGRA parameter 
space, and thus our quantitative results should not be taken too precisely. Choosing the central 
values of the neutrino parameters keeps the presentation simple. With these assumptions, the 
diagonalizing matrix becomes 

( .83 .56 \ 
-.4 .59 .71 
.4 -.59 .71 



U, 



\ 



(A.4) 



I 



For the normal hierarchy we have m Ul « m„ 2 <C m U3 = M, while for the inverted hierarchy we 



have m U3 <C m Ul r 
We then find 

r(x" 



m U2 = M. 



M 2 \U i3 U* 3 \ 2 



Normal Hierarchy 



M 2 \Un U* x + U i2 U* 2 1 2 Inverted Hierarchy 



(A.5) 



The branching ratio BR(x + — > (\vj) = T{x + ~^ ^t u j)/^(x + ~^ an y is given by mf- / J2ki m 
which is easily evaluated as 

( o \ 

25% 25% 
25% 25% 



2 



BR( X + -> g+ Vj ) 



V 



for the normal hierarchy, 



(A.6) 



I 



BR(x" 



4< 



for the inverted hierarchy. 



(A.7) 



( 50% o N 

12.5% 12.5% 
\ 12.5% 12.5% J 

In the case of final state neutrinos, one sums over all flavors since they are not detected. All of the 
above holds for the analogous decays involving SUSY partners: \ + ~^ ^t v 3 an< ^ a ^ so wm Si ye us 
the flavor structure when, for example, a vr of flavor j decays into a x + * an d a lepton t~ 



References 



[1] For a review with an extensive list of references, see M. Sher, Phys. Rept. 179, 273 (1989) and J.F. Gu- 
nion, H.E. Haber, G.L. Kane, and S. Dawson, The Higgs Hunter's Guide (Westview Press, Boulder, 
Colorado, 2000). 



17 



[2] H. E. Haber and G. L. Kane, Phys. Rept. 117, 75 (1985). 
[3] T. P. Cheng and M. Sher, Phys. Rev. D 35, 3484 (1987). 

[4] R. M. Barnctt, G. Senjanovic, L. Wolfenstein and D. Wyler, Phys. Lett. B 136, 191 (1984); R. M. Bar- 
nett, G. Senjanovic and D. Wyler, Phys. Rev. D 30, 1529 (1984); Y. Grossman, Nucl. Phys. B 426, 
355 (1994) |arXiv:hep-ph/9401311| ; V. Barger, H. E. Logan and G. Shaughnessy, Phys. Rev. D 79, 
115018 (2009) [arXiv:0902.0T70l [hep-ph]]. 

[5] H. S. Goh, L. J. Hall and P. Kumar, JHEP 0905, 097 (2009) |arXiv:0902.081~4l [hep-ph]]; M. Aoki, 
S. Kanemura, K. Tsumura and K. Yagyu, arXiv:0902.4665 [hep-ph]; S. Su and B. Thomas, Phys. Rev. 
D 79, 095014 (2009) |arXiv:0903.0667l [hep-ph]]. 

[6] S. Gabriel and S. Nandi, Phys. Lett. B 655, 141 (2007). 

[7] S. M. Davidson and H. E. Logan. larXiv:0906.3335l [hep-ph] . 

[8] C. D. Froggatt and H. B. Nielsen, Nucl. Phys. B 147, 277 (1979). 

[9] A. Aranda and M. Sher, Phys. Rev. D 62, 092002 (2000) |arXiv:hep-ph/0005113] . 
[10] M. Spira, A. Djouadi, D. Graudenz and P. M. Zerwas, Nucl. Phys. B 453, 17 (1995) 

[arXiv:hep-ph/9504378] . 
[11] J. Alwall et al, JHEP 0709, 028 (2007) [arXiv:0706.2334 [hep-ph]]. 

[12] J. Pumplin, A. Belyaev, J. Huston, D. Stump and W. K. Tung, JHEP 0602, 032 (2006) 

[arXiv:hep-ph/0512167] . 
[13] G. Aad et at [The ATLAS Collaboration], larXiv:0901.0512l [hep-ex], 

[14] A. Freitas, A. von Manteuffel and P. M. Zerwas, Eur. Phys. J. C 34, 487 (2004) |arXiv:hep-ph/0310182] . 
[15] We are assuming here that its supersymmetric partners are not much lighter-if so, other two-body 
decay channels which are observable could open up 



18 



